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EIGENVALUE ESTIMATES FOR A THREE-DIMENSIONAL 
MAGNETIC SCHRODINGER OPERATOR 
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O 

J 1 Abstract. We consider a magnetic Schrodinger operator i/" = (— i/iV — A)^ 

frt ' with the Dirichlct boundary conditions in an open set O C K^, where h > is 

a small parameter. We suppose that the minimal value bg of the module \B\ 
of the vector magnetic field B is strictly positive, and there exists a unique 
^y. • minimum point of \B\, which is non-degenerate. The main result of the paper 

' is upper estimates for the low-lying eigenvalues of the operator H^ in the 

semiclassical limit. We also prove the existence of an arbitrary large number 
of spectral gaps in the semiclassical limit in the corresponding periodic setting. 

"^ , 1. Preliminaries and main results 

1.1. Main assumptions. We would like to analyze the asymptotic behavior, in 
the semiclassical regime, of the low-lying eigenvalues of the Dirichlet realization of 
the magnetic Schrodinger operator in an open set O in K^: 

H'^ = {hDx, - Ai{X)f + (hDx, - A2{X)f + {hDx, - A^{X)f , 

CN \ where A = (Ai, yl2, A3) e C°°(f2,R'^) is a magnetic potential and /i > is a small 

^-p ■ parameter. We will denote the coordinates in R^ as X = {Xi,X2,X^) = {x,y,z) . 

Let B = rot A = (i?i, i?2, -B3) be the corresponding vector magnetic field: 

O' B =^-^ B = — ~^ B ^—~ — 

CN ' dy dz ' dz dx ' dx dy 

Put 

60 =mm{\BiX)\ : X en}. 

^^ I We assume that there exist a (connected) bounded domain fii CC il and a constant 



H ■ £0 > such that 



(1.1) \BiX)\>bo + eo, x^Oi. 

We also assume that: 

(1-2) 60 >0, 

and that there exists a unique minimum Xq G 51 such that |i?(Xo)| ~ bo, which is 
non-degenerate: in some neighborhood of Xq 

(1.3) C-^\X - Xof < \BiX)\ - 60 < C\X - Xof . 
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1.2. Main statement. 

Denote 



d = detHess|B|(Xo), a= —^{}iess\B\B ■ B){Xo) . 
Denote by Xo{H'') < Ai(iJ'') < X2{H'') < ... the eigenvalues of the operator H'^. 



Theorem 1.1. Under current assumptions, for any natural m, there exist Cm > 
and hm > such that, for any h G (0, hn 

1/2 

(2m + l) + v 



\m{H'') < hbo + h^l'^a^l'^ + h^ 



1 / d 

'^ V 2^ 



+ c^h^/' 



where v is some explicit constant (which means that it is given by a rather compli- 
cated explicit formula). 



Theorem 1 1 . II wiU be based on a construction of quasimodes. 

Theorem 1.2. Under current assumptions, for any natural j , k andm, there exist 
'^j,k,m ^ (^'i^i^)^ Cj^k,m > and hj^k,m > such that 



(1.4) 


i^'nM 


"il ' ^J2,k2,m 


J 


and, for 


any h G 


(0,/lj,fe,rn], 




(1.5) 




II Tjhih 


- 


where 









^jlJ2'^kik20mim2 



o 



f'j,k,nAj,k,. 



<c. 



j.k.jn 



ji ^ki ,m-i ,J2 ,k2 ,m2 



hHct^lk,J 



(h). 



with 



f-j,k,m — f-j,k,m,Qh + /^j,fc,m,2^^ + f'jMjn,4h . 

H,k,m,Q = {2k + l)6o , l^j,k,m,2 = (2j + 1)(2A: + ly/^a'/^ , 



1/2 



fJ'j,k,rn,4 

where v{j, k) has the form 

v{i,k) = V22{2k + lf 

with some explicit constants I'd, i^ii, ^^22- 

Remark 1. 

It is conjectured that 

■" 1 



^(^^) (2m+l)(2fc+l) + z.(j,fc). 



J^ii(2j + 1) 



vg 



\n{H'') >hbo + h^'^a^'^ + h^ 



2bn \2a 



1/2 



{2m + l) + iy 



Crah^" 



1.3. History of the problem. 

May be at the level of the mathematical problem, the starting reference for the 
spectral analysis of self-adjoint realizations of the magnetic Schrodinger operator is 
the paper by Avron-Herbst-Simon [T] where the role of the module of the magnetic 
field in the three-dimensional case appears for the first time. Further investigations 
were inspired by R. Montgomery [T3], who was asking "Can we hear the locus of 
the magnetic field" (by analogy with the celebrated question by M. Kac). In [Hj . 
this question was studied for the two-dimensional magnetic Schrodinger operator. 
Motivated by the question of R. Montgomery, the first author and Mohamcd in [12] 
investigated the asymptotic behavior of the low-lying eigenvalues of the Dirichlet 
realization of the magnetic Schrodinger operator in the case when the magnetic 
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field vanishes. This study was continued more recently in [TH [71 [H |3] (see also 
[8]). The case when the magnetic field never vanishes was analyzed in detail for the 
Dirichlet realization in the two-dimensional case in [TJ] and more recently in |10[|llj . 
Moreover, there is a big literature devoted to the spectral analysis of the Neumann 
realization because of its connection with problems in superconductivity (see [4] and 
the references therein). Finally, we do not also give a complete description of the 
semi-classical results obtained in the case when an electric potential V is creating 
the main localization and refer to [S] and [5] for a presentation and references 
therein. 

The paper is organized as follows. In Section [2l we make some simplifications of 
the operator H'^, using linear changes of variables and gauge transformations, and 
present the normal form appearing in the generic situation. Section [3] is devoted 
to the analysis of the action of the metaplectic group on our models and to the 
introduction of a suitable rescaling of the problem. In Section 2] wc construct ap- 
proximate cigenfunctions of the operator H^, completing the proof of Thcorcm ll.2l 
In Section [5] wc consider the case when the magnetic field is periodic. Wc combine 
the constructions of approximate cigenfunctions given in Sections |4] with the results 
of [6] to prove the existence of arbitrary large number of gaps in the spectrum of 
the periodic operator 7J'' in the semiclassical limit. 

2. Preliminaries for the proof: Normal form 
In this section, we make some simplifications of the operator _ff'', using linear 
changes of variables and gauge transformations. 

Without loss of generality we will assume that Xq = (0,0,0) and the magnetic 
field B at Xq is (0, 0, bg) . Thus, we can write the Taylor expansion up to order 4 
at (0, 0, 0) of the magnetic field: 

(Bi{X)=h-X + Qi{X) + Ci{X) + Ri{X) + Oi\Xf) , 

(2.1) I B2iX) =e2-X + Q2{X) + C2{X) + R2{X) + 0{\X\'>) , 

i BsiX) = 6o + QziX) + C3(X) + i?3(X) + 0{\X\^) , 

where ij € S.'^,j ~ 1, 2, and Qj, Cj and Rj, j = 1, 2, 3, are the terms of order 2, 3 
and 4, respectively. 



By assumption (|1.3p . we have 

(2.2) {£i ■ Xf + {£2 ■ Xf + 2boQ3{X) > , VX e 



We also have 



which implies 



, - dBi dB2 dBj, 
ox ay oz 

d{£i-X) ^ d(i2-X 



(2.3) 



dx dy 

dQi ^ dQ2 ^ dQs 

dx dy dz 

dCi dC2 dC3 

dx dy dz 

dRi ^ dR2 ^ dR3 

dx dy dz 



= 0; 
= 0; 
= 0; 
= 0. 
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With the notation, 

^1 = (ai, 02,^3) , -^2 = (/3i,/32,/33), 
the first hne of p.3p gives 

(2.4) ai+/32=0. 
We can assume without lost of generality that 

(2.5) a3=0. 

Indeed, either it is already the case, either /33 = and we can exchange the role of 
X and y, or a^ ^ and (3^ =/= and a rotation around the z-axis permits to obtain 
the cancellation of 03 . 
With the notation 



(2.6) Qi{X) = Y, a^JX^XJ , Q2{X) = ^ h,X,Xj , Q:i{X) = ^ c^X^X 
j,j=i »,i=i ij=i 



3 3 

' ^. . V. Y . 



the second line of (|2.3p reads 

(2.7) aij+52j+C3j =0, .7 = 1,2,3. 
With this notation, (jl.3p is equivalent to saying that the matrix 

(2.8) = iHcss(|Bn(0) = 6oHess|B|(0) 
or equivalently (cf. (|2.2p ) 

/ af + /?? + 2&0C11 aia2 + /3i/32+26oCi2 /3i/33 + 26oCi3\ 
Q = aia2 + /3i/32 + 260C21 «! + /3| + 260C22 P2Pz + 260C23 
\ /3i/33 + 260C31 ;S2/S3 + 260C32 Pi + 260C33 / 

is positive definite. We will denote the elements of Q by {qij}i,j=i,2,z- 

Without loss of generality, we can, after possibly a gauge transformation, assume 
that: 

Ai = 0. 
With this assumption, we get an equation for A2 by writing: 

S3 = ^ = &0 + QsiX) + C3(X) + i?3(^) + O(l^l') , 

and without loss of generality, we can, after possibly a gauge transformation assume 
that 

^2(0,y,z) = 0. 
This implies 

(2.9) A2{x,y,z) 

pX pX pX 

^box+ Q3iX)dx+ C3{X)dx+ RsiX) dx + xO{\X\^) , 
Jo Jo Jo 

where X = (x, y, z) . Similarly, we get an equation for ^3 by writing; 

B2{x, y, z) = —^ ^i2-X + Q2{X) + C2{X) + R2{X) + 0{\X\^) 
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and without loss of generality we can, after possibly a gauge transformation (keeping 
the previous properties) assume that 

A3(0,0,z) = 0. 

This implies 

nx px px 

A3{x,y,z)^~ h-Xdx- Q2{X)dx- C2{X)dx 

(2.10) -^^ -^o ■'' 

- R2{X)dx + a3{y,z) + xO{\X\''), 
Jo 

where X — (x, y, z) and 03(1/, z) — ^3(0, y, z) satisfies 

(2.11) a3(0,z) = 0. 

From (|2.9p and (|2.10p , we obtain, using the definition of Bi , 

B,ix,y,z)=^(X)-^iX) 
oy oz 

= -h^~l — (X)dx-y^ -^(X)dx^l —iX)dx 

no 

^{X)dx + xO{\Xf), 

and using (j2.3p . we obtain 

Bi(x, y, z) =aix + Qi{X) - gi(0, y, z) + Ci(X) - Ci(0, y, z) 

+ R^iX) - i?i(0, y, z) + ?^(y, z) + 0(|X|^) . 

oy 

On the other hand, comparing with (j2.ip . we obtain 

^(y,z)=a2y + Qi(0,j/,z) + Ci(0,j/,z) + i?i(0,y,z) + O(|y|5 + |z|5). 
Thus, we obtain, having (j2.1ip in mind, 

(2.12) a3(y,z) = ia22/' + y" Qi{0,y,z)dy 

+ f Ci(0,y,z)dy+ f'\iiO,y,z)dy + yO{\yf + \zf). 
Jo Jo 

Coming back to (|2.9p . we have 
^2(2;, y, z) =&o2; + -ciix^ + ciax^y + Ci32:2z + C22xy^ + 2c2,xyz + C33xz^ 

pX pX 

+ / C3(X)dx+ / R:i{X)dx + xO{\X\''). 
Jo Jo 

Similarly, coming back to (|2.10p and using (|2.12p . we obtain 

^3(2;, y,z)^- 2^^^^ ~ ^^^^ ~ ^^^^ '^ 2"^^^ 

+ / Ci(0,y,z)dy+ / Ri{0,y,z)dy 
Jo Jo 
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- ^&iia;^ - 612x^2; - bi^x'^z - 6222:2/^ - '2b23xyz - bs^xz^ 

+ :;a22y^ + a23y'^z + assyz'^ - / C2{X)dx- / R2{X)dx 
■J Jo Jo 

+ xO{\X\') + yO{\y\' + \zf). 
Finally, what we have got is the following normal form: 

Proposition 2.1. Using only linear change of variables and gauge transformations, 
we can assume that 

Ai{x,y,z)=0, A2i0,y,z) = 0, A3(0,0,z) = 0. 

Moreover: 

A2{x, y, z) ^box + -ciix^ + c^x'^y + ci^x'^z + C22xy^ + 2c23xyz + c^^xz^ 



/ C3{X)dx+ / R3{X)dx + xO{\X\'^), 
Jo Jo 



^3(2;, y,z) =~ -fiix'^ - I32xy ~ jS^xz 

1 py py 

(2-^^) +:^a2j/'+/ C,{0,y,z)dy+ R,{0,y,z)dy 

^ Jo Jo 

- -^wx^ - bi2x'^y - bi3X^z - b22xy^ - 2b23xyz - 633x2^ 

+ ::a22y^ + a23y^z + a33yz^- / C2{X) dx - / R2{X)dx 
-J Jo Jo 

+ xO{\X\^) + yO{\yf + \zf). 

3. Preliminaries for the proof: Metaplectic transformations 

In this section, we make further simplifications of the operator H^' , using meta- 
plectic transformations. First, we introduce a suitable ft,-dcpendcnt rescaling of the 
problem and expand the resulting operator in fractional powers of h. 
3.1. Rescaling. 
Let us make a change of variables 

x = h^x, y = h^y, z = h^z. 

We will only apply our operator H^ on functions which are a product of cut-off 
functions localized in a neighborhood of (0,0,0) £ R"^ with linear combinations of 
terms like 

h''w{h-'/^x,h-^^^y,h-^^^z), 

with w in 5(R'^). These functions are consequently 0{h°°) outside a fixed neigh- 
borhood of (0,0,0). We will start by doing the computations formally in the sense 
that everything is determined modulo 0{h°°), and any smooth function will be 
replaced by its Taylor's expansion at (0, 0, 0) . We introduce 

Ci{x) ~ 2_^ PijkXiXjXk , 

i<j<k 

C2ix) = 2_^ qijkXiXjXk , 



EIGENVALUE ESTIMATES FOR A THREE-DIMENSIONAL MAGNETIC OPERATOR 



i<j<k 

We will also need the coefficient of z^ in Rj : 

Rj{X)^S,z^ + ..., j = l,2,3. 



We have 



where 



H'' = hP^ , 



P'' = Di 



(Oy - h-^/^A2ihix, hiy, his)) 

Now, using (j2.13p . we expand the operator in fractional powers of h. First, we 
compute 

= b^x + h^cssxz'^ + hi [cisx'^z + 2c23xyz + rss^xz^] 

+ h -ciix^ + ci2X^y + C22xy^ + -risax'^z^ + r2s.zxyz^ + 5s.x 

+ 0{h^'^) , 
and 

h^^/'^A3{h^x, h^y, h^z) 

= - hi P^xz + h^ 



2^^ 



1 



r.z2 



hiy + 2'^22/ - h3.S:z + aggyz 



+ hi [-bi^x^z - 2b23xyz + a23y^z - qasaxz^ +P333yz^] + 0{h) . 
Using these formulae and omitting the tilda's, we obtain (in the dilated coordinates) 

P^ = P(, + h-iPi + h^-P2 + hiP3 + hPi + 0{h^l^) , 
where 
Po =Dl + (Dy - hoxf , 
Pi=0, 

P2=- 2c33XZ^{Dy ~ box) + {D, + P3xzf , 

P3 = - 2ci3x'^z{Dy - 602;) - 2c23[{Dy - box)xyz + xyz{Dy - b^x)] 
- 2r333Xz^{Dy - bax) + {Pix^ + 2hxy - a2y^){D, + f33xz) 
+ {b33X - a33y){{D, + p3xz)z^ + z'^{D, + /^gxz)) , 



P4 = - {Dy - box) 



-cnx^ + ci2X^y + C22xy^ + -r\33X^z^ 



T233Xyz 



S3XZ 



-ciix^ + ci2X^y + C22xy'^ + -ri33X^z^ + rzasxyz^ + 63xz^ 



{Dy - box) 



+ clr^x'^Z^ + {D^ + I33XZ) [6132:^2 + 2b23Xyz - a23y^Z + q333XZ^ - P333yz^] 

+ [bi3X z + 2b23xyz - 0231/ z + 9333x2 - P333yz \ {D^ + P3XZ) 
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3.2. Partial Fourier transform and gauge transform. 

Now wc arc going to apply some nictaplcctic transformations. Wc recall that these 
transformations are unitary and therefore preserve the spectrum of the operators. 
First, wc make a partial Fourier transform in y: 

Pq =Dl + (r/ - b^xf , 
Pi=0, 



P2 

Ps 



Pi 



- 2c33xz^{r] - box) + (D^ + Psxz)^ , 

- 2ci3X^zir] - box) - 2c23[(?7 - bQx)x{-D,,)z + x{-D^)z{ri - b^x)] 

- 2r333Xz^{T] - box) + {/3ix^ + 2p2x{-D^) - a2D^^){D, + ji^xz) 
+ (633a; + a33A,)((-D^ + P-,xz)z'^ + z'^{D^ + P^xz)) , 



rl 



{rj - box) ^CliX^ + Ci2X^{-Dr,) + 022X0^^ + -ri33X^2^ 



+ r2332^(-A,)z +5:ixz 



-ciix '\'Ci2X {-D,^)+C22xD^ + -ri3,3 



X z 



+ r2'i?,x{--Dn)z^ + S^xz"^ (77 — box) + 033X^2^^ 

+ {Dz + I3?,xz) \bi'ix^z + 2b23x{-Djj)z - a23D'^z + qsssxz^ - Psdsi-DTj)^^] 
+ [bi^x'^z + 2b23x{-Dn)z - a23D^z + qzs.3xz^ - P33.3{~Dy^)z^] {D^ + Psxz) 

+ ( -Pix^ + hxi-Dr,) - -a2Dl + b^s.xz'^ + a33£'„2^ j • 

Next, we make a translation x = x — ■?- (and forget the "check"): 



2^2 



Po =Di + b'ox 
Pi=0, 

P2 =2c33(x + ^)z'-box + {D, + I3i{x + ^)z)2 , 
bo bo 

P3 =2C13(.T + ^)^z{box) + 2C23 






D^- Dnjz 



+ (^+ T-)(-r^^ - Dv)4box)] + 2r333(x + ^)z^{box) 
bo \bo /J 60 

+ hi{x+^)^ + 2f32{x + ^)(^D, 



D, 



a2 1 ■j-P'x - Dr, 



X (i?, + /33(X + ■^)z) + (633(2; + J-)- 033 (7- A - A, 
Oo ^ Oo V60 

X ((A + /33(a; + ^)z)z^ + ^^(A + /33(a; + j-)z)) , 
V 60 00/ 

-~box\-cii{x + ^)^ + ci2(a- + -^)2 ( —A - A, 

Lo Oo Oo \Oo 



+ C22(a: + ^) (—A- A 
Oo Voo 



+ ^'^133(2: +;J^)'22 
^ Oq 
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A Oq Oo \Oo 

+ ^233(2; + r-) f 7^^- - ^')) ^' + -^alx + ^)^^] (feoa:^) 



60^ 



+ 4(^ + f )'^* + (^. + /33(x + ^)z) 6i3(x + f )^z 



n 

bo' 



bo' 



+ 2623(2; + y) ( ^Dx -DAz~a23i yD^ - A, I z 
+ 9333(2; + ^)z^ - P333 lyDx - dJ z3 



'?^2. 



^nH 



60' 



^0 V^o 



^13(2; + r^Yz + 2623(2: + 7^) -D.X -DAz- 023 ^D, - D 



1 



bo 



ZL^~3 

60' 



1 

b'o 



9333(2; + f)z^-p333( fi?, -A, )z' (A + /33(2; + f )2) 



60' 



(^^^(^• + ^)^+^^(^ + ^) 



60^ 



feo 



-'2; J-^ri 



'' ^ ' ^-A - Aj - ^^2 I ^A - A 

+ &33(2;+^)z^-a33 (— A- A )Z^' 

Oo \oo J 



Finally, we make a gauge transformation (by exp(— i^Tyz^)): 



Po =Dl + 6^x2 , 
Pi -0, 



'^.2 



P2 =26oC332;(a; + -!-)z' + (A + iS^xz) 
bo 



^ \2 



V 



P3 =2boCi3xix + -^Yz + 2c23 (boxYx + f ) -A + ^z' -DAz 



bo 



bo V60 



/3: 



3 .2 



26o 



A 



'3 -,2 



^^-,3/ 



(•^ + I^) 7-^- + :^^ - Dr, z{hox) + 2r333(2; + ■^)z'\box) 



bo'ybo"-' ' 26o~ " 



bo' 



hiix + ^Y + 2P2ix+^) 
V 60 60 



I 1 n , ^3 2 n 

"2 I 7- A + 777-z - A 

Oo ^Oo 



—A + —z^ - A 
60 26o " 



(A + ^32:z) 



'/ X „ / 1 n , /^3 9 



+ &33(2^ + r) - «33 T-l?x + ^z' - Dr, 

' 60 V60 26o 



X ((A + P3xz)z'^ + z^{D^ + I33XZ)) 
Pi =box 



1 (x+f)-Vci2(x+f)^(lA , „, 

3 60 60 Voo 26o 



^^.^-A 
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2 






+ r233(x + ^) ( ^i?. + ^Z^ - A, ) Z^ + <^3(x- + ^)z^ 






rl ' '^3 , „ /^, , ^/n2 n r. , /^3 2 



6o 

r233(:^ + r) (r^- + 1^^' - ^') ^' + '^3(X + ^)z4l (box) 

bo V5o 26o / oo J 



'?^2 4 



'?n2. 



+ C33(a; + 7^)'z^ + (A + I33XZ) bi3{x + -^Yz 



ho' 



'? \ / 1 



A 



'3 „2 



2623(x + 7^) T-Dx + T^z^ -DAz~ a23 j-D, + ^z' - D,, z 
bo \bo 2bo J \bo 2bo 



^3 2 



'?^,3 



bo' 



/3: 



3 _2 



9333(3^ + 7^)^ -P333^A + :^Z -A, U^ + 6i3(x + f)^Z 



bo 2bo 



'?^2. 



&o' 



A^ 



3 „2 



2b23{x + 7^) -A + T^Z^ - A M - ^23 r A + ^z' " A 

k .2 



/3: 



3 ^2 



+ 9333(2;+ —)z^ -P333 ( 7--Dr + ^Z-' - Dr^ I Z 
Oo V bo 26o 



(i^, + (33xz) 



(biix+^r+i32ix+^) 

V 2 60 Oo 



^ N /^ 1 n , ^3 2 n 
60 26o ' 





1 ^1 " /^3 2 r^\\ , f^ , V.2 



- 033 (^D. + ^z^ - A ) z^ 
\ Oo ^oo I 



4. Construction of approximate eigenfunctions 

In this section, we complete the proof of Theorem 11.21 We start with a formal 
construction. We will construct an approximate eigenfunction of the operator P'' 
in the form 

corresponding to an approximate eigenvalue 

We will express the cancelation of the coefficients of h!'!^ {£ £ N) in the formal 

expansion of (P'' — A'')?i'', starting with 1 = 0. 

4.1. The first equations. 

The first equations read (here we omit tilda's): 

(4.1) (Po - Ao)iio =0 , 

(4.2) (Po-Ao)ui=Aiuo, 
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(4.3) (Po - Ao)u2 = - P2U0 + Aiui + A2U0 , 

(4.4) (Po - Ao)m3 = - P3U0 - P2U1 + X1U2 + A2M1 + A3M0 , 

(4.5) (Po - Ao)u4 = - P4U0 - ^3^1 - P2U2 + A1U3 + A2U2 + A3U1 + A4M0 , 

and we now show how to solve them successively. 
4.2. Main equation. 

The first equation (|4.ip reads 

{Dl + blx" - Ao)uo = . 
We arrive at the eigenvalue problem for the harmonic oscillator 

(4.6) [,1 := Dl + by . 
Recall that its eigenvalues have the form 

Sp([)i) = {/ifc = (2fc + l)6o:fceN}. 
An eigenfunction of t)i associated with the eigenvalue ^k is given by 

where H^ is the Hermite polynomial: 

i^.(0 = (-l)V^^(e-^). 

The norm of (l)k in L^(IR, dx) equals the norm of H^ in i^(R, e~* dt) , which is 
given by 



Some well-known formulae, concerning to the Hermite functions, which we need in 
the paper, are gathered in the appendix. 

We will take Aq as the fc-th eigenvalue of ()i 

(4.7) Ao = (2fc + 1)60 
and will look for a solution uq to (|4.ip in the form 

(4.8) uo{x,T],z) = (l)k{x)xoiv)'^n{z), 

with xo and ^Pq in 5(R) to be determined in the next steps. 

Remark that what we do at this step corresponds to the so-called reduction to 
the fc-th Landau-Level. 

4.3. Second equation : coefficient of ft,4. 
We look for a solution ui to (14. 2p in the form 



(4.9) ui{x,'q,z) = (j)k{x)vi{r],z), 

with vi G 5(M^) to be determined in the next steps, and take 

(4.10) Ai=0. 
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4.4. Third equation: coefficient of /i^. 

In order to solve the equation (j4.3p . we find as a necessary condition that, for any 
z we should have, 

(4.11) - ^:^(-P2Uo, <^fe). + A2Xo(ry)*o(^) - , 

with (•, ■)x denoting the scalar product in l?{^x). Using that 4)\ is even, we obtain 
that 

(4.12) Wk\V WkV 

=Xo(?/)f)3*o(^), 

where 

2fc + 1 

(4.13) [)3:-i?,' + ^r^(/3l + 2c336o)z2. 

Thus, the condition (|4.1ip can be rewritten as 

(4.14) fl3*o = A2*o. 

This determines A2 and ^0. We take \i as the j-th eigenvalue eigenvalue of (53: 



(4.15) A2 = (2j + l)(2fc + 1)1/2 J^i±^ ^ (2j + 1)A2 

where 

/g33 



(2fc + 1)1/2 S^2^=(2fc + 1)1/2^ 
V 2oo 



and 'I'o as the corresponding eigenfunction: 

M/o(z) - ^,(z) := A^/'i/,(A^/^z)e-^^^V2 , 
Observe that 



933-:^^^(Xo)-25oa. 



192|B|2 

2 az2 

Thus we obtain 

A2 = (2j + l)(2A: + l)i/V/2. 

To find U2, from (|4.3p we obtain that 

([)i - Ao)u2 

= - (l>k{x)xo{v)Dlipjiz) - a;0fe(a;)xo(77)(2c33'7z^ + ^DzZ + zDz))ipj{z) 

- {Pa +'2boC33)x'^(l)k(x)xo{v)z^ip3{z) + ^2<l>k{x)xo{v)i'3{^) 

= 775 i(l)k+i{x) + 2k(f>k-i{x)) xo{v){2c33Vz^ + h{D,z + zD,))i;j{z) 

2V 

- -^[(il + 260C33) (0fc+2(a;) + 4fc(/c - l)(/)fc-2(a:)) xo(??)^'V'j-(z) . 

Therefore, using the identity 

(f]i - Ao)0fc+f = 2ibo4'k+e, 
we get 
(4.16) U2ix,r/,z) = U2{x,r],z) + cj)k{x)v2{ri , z) , 
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where 

U2{x,r],z) 



46; 



3/2 



bk+i{x) - 2fc0fc_i(x)) xo(f?)(2c3377z2 + PaiD.z + zD^))i>j{z) 



ml 



{Pl + 260C33) {.^k+2{x) - Ak{k - I)4>k-2{X)) X0{V)Z^4'J{Z) : 



and V2 e 5(K^) has to be determined in the next steps. 
4.5. Fourth equation: coefficient of /i*. 

In order to solve the equation (j4.4p . we find as a necessary condition that, for any 
{z,ri) we should have 



(4.17) 



1 



\Hk\ 



(-P3U0, 



1 



'Jk/x 



\Hk\ 



{P2Ui,(t)k)x + >^2Viiz,r]) + X3Xoi'n)'^j{z) =0. 



For the first term in the left-hand side of (|4.17p . since (pl is even, we have 

{P3Uo,(l>k)x = ■7-TT-r^{P:^U0,(l)k)x, 



where P^ is the even part of P3 as a differential operator in a;. As a sum of 
homogeneous components in (77,!)^), the operator P^ is written as 

where 



^i'^ = (^1^ - 2/^2^ A, - a2Dl) D,, 



P 



(1) 



„2 '/ , „ I ^7 /„n , n ^\ o^2r- '^ " ^ 



29132; 1- + 923 -pii^^x + ^x^) - 2a; D^ + a2l33 — {xDx + Dxx)D^ 



^2^3 , "2/^3 r, , u '^ , n 



2b. 







2&I 



'0 



50 



{z^D, + D,z'), 



3(0) _ 923/33 2^3 



p^,u; ^^^^2^3 ^ 2(/?3fe33 + 6or333)x2z3 + {^^x^ + 2/32-^xD, - a2^Dl)D, 

Oo Oo Og 



"2/31 , /33a33 



2&n 



(D,z4 + z^D,) 



"2/?! , /33a33 



26g 



[xD^ + D^x)z^ 



(4.18) 



Now we compute: 
1 



\Hk 



{P3UQ,(j)k)x 



1 



where 

T2(^3^^"O,0fe)2^ 



3(2) 



pr^o, 



\HkP''^ 



.(1) 



K^^"o, 



P|°''wO,(/'fc)x 



l^^l 



/3i|j - 2/32^A, - «2A' ) Xo(^)i?.^j(2) 



("/Si^ - 2/?2^A, - a2i?^) Xo(r/)iA^/'z (^j+i(^) - 2j^j-i{z)) , 
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JJ^iPt^^o.Ah 



= (2fc + l) 



V 



1 



91372 ~ <l2zi-D^ xo(»7)z'0j(z) 
^2/33 , ^2/33 n 1 ?, ^ , n 



26 







= (2A: + 1 



26o 

J] 1 

■ — /7o ■:> — 

^0 

132 k , a2/33 



■60 



Xo(77)(2'^.+^.^')'/'jW 



913 72 ^ 923 — ^1; 



1 



Xo(?7)— "ITI W'i+i(^) + 2jV'i-i(2)) 



2A 



26g 



77 1 1 

■V + -ITT—Dr, + 6337- + a33-D^ Xo('7) — rrt*['0i+3(z)+ 
^oo oo J 4A2' 



+ (2j + 2)^,+i(z) - 4/V;,_i(z) - 8j(j - l)(j - 2)V^j_3(^)] , 



^ (Pi°^uo,0fe). 



IfffcF' ' 



= (2fc + l) 



923/33 , /33633 + 60^333 



262 



60 



Xo(?7)z^V'j(2) 



,^ 2fc + l ^ 1 . 2fc + l, . , , , , 

(Pi^^T — + P2 7-« - Q^2 „, )xa{'n)Dz4'Az) 

ra2/3| , /33a33 



Xo(r;)(z4i?,+i?,z4)^,(z) 



= (2fc + l) 



26o 

923/33 , /33633 + ^0?'333 



26g 



Xo(?7) 



8A^/^ 



(^5+3(2)+ 



+ (6j + 6)7/^,+! (z) + 12j^V^,-i(z) + 8jij - l)(j - 2)7A,_3(z)) 

^ + ^1 XoW^^*(V',+5(^) + (6j + 12)V',+3(^) 
8&0 26o J I6A2' ^ 

+ 4(2j2 + 4j + 3)^j+i(2) - 8(2j2 + l)ji^,_i(z) 

- 48j(j - 1)'(.? - 2)^,_3(z) - 32j(j - 1)0- - 2)(j - 3)(j - 4)^,_5(^)) • 

For the second term in the left-hand side of (|4.17p . as in (|4.12p . we obtain 

1 



\H,P 



{P2Ui,4>k)x = ^3Vi(ri,z) 



Thus, the condition (j4.17p can be written as an equation with respect to vi in 
the form 

1 



(4.19) 



(f)3 - X2)vi{ri,z) 



\H, 



{P^uo,<j)k)x + X3Xo{v)'fPjiz) 



In order to solve (|4.19p . as a necessary condition, we should have 
1 



iHkPWH, 



r(-f3 ""0, (t>k^])x,z + AsXcC??) = : 
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with {■,-)x,z denoting the scalar product in L^(]R^_j). It follows from (|4.18p that 

s odd in 2 

A3-O. 



the function {P^uq, 4>k)x{V! ^)4'j{^) is odd in z, therefore, {P^uq, 4'k4'j)x,z ~ and 



(4.20) 

Using (|4.18|) and the identity 

(()3 - X2)ipj+e = 2iA2'ipj+e, 
we find a solution to (|4.19p in the form 

(4.21) 
where 



(2) , (1) , (0) 



,(2) 



,(1) 



f/?i^ - 2/32^A, - a2D^) ;^o(^)^^i {^j+i{z) + 2j^,_i(z)) , 



913 72 -923^ A, 

On OO 



2fc + 1 

Xo(??) 3/2 (V'i+l(^) - 2jV^i-l(2:)) 



4A, 



— —^ + -;^7^ A + 033— + 033 A 



26^ 



26o 



•feo 



Xo(?7) 



24A 



3/2 



i[V'j+3(^;)+ 



+ (6j + 6)V,+i(z) + 12jV,-i(2) + 8j(j - l)(i - 2)^,-3(2)] , 



,(0) 



'?23/33 , /33fe33 + &0?'333 



2&2 



fco 



Xo('7)-r77H(^J+3(^) + 



48A 



5/2' 



+ (18j + 18)V;j+i(z) - m'fi>,-M ~ 8jij ~ l)(j - 2)Vj-3(^)) 

2fc + 1 1 2fc + 1 1 

~ (/^i oA + /^2 7-i - "2^77 — )Xo(f])— n7^i(V'j+i(^) + 2jijj-iiz)) 
zoq Do Zbo 4A, 



"2/31 , ^3033 



+ 



2&0 



Xo(^) _.5/2 ^(y^V'3+5(^) + (.?■ + 2)l/;j+3(z)- 



16A 



+ 2(2j2 + 4j + 3)^,+i(z) + 4(2j2 + l)j^,_i(z) 
+ 8j(j-l)'(j-2)^,-3(^) 



2/. o^.,. /.^ : 16j(j--l)(j-2)(j-3)(j-4) 



V'j-5(^) 



With such a choice of wi, the orthogonality condition (|4.17p holds, and a solution 
U3 to (|44)) exists. Substituting (|4.2ip in (|49l) . we find wi. 
4.6. Fifth equation: coefficient of h. 
The necessary condition for solvability of (|4.5p reads 



(4.22) 



U?fc 



r{P4:Uo,4'k)x 



Hki 



;{P3Ui,(l)k)x 



1 



(P2U2, 



Using (|m|) and (|i?T^ . we get 
1 



^ |j2(i^2W2,(?!'fc)x 



+ ^2(^2, 0^)2; + A4Xo('7)'0j(^) = 0. 

:{P2U2,(t>k)x + ()3W2 (77,-2), 



and 



\Hk¥ 



\Hk 



{u2-,4>k)x = V2{ri,z) . 



Therefore, the condition ()4.22p can be rewritten as an equation with respect to V2 
in the form 
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(4.23) H}3-X2)v2iv,z) 

\\-tik\\ \\-"k\\ ll-"felr 

The necessary condition for solvability of (|4.23p is given by 

{P2U2, (t)ki>j)x,z + A4Xo('7) = . 



Lemma 4.1. The identity (|4.24p has the form 

I A o^\ 2^ + ^ "^22933 - 9^3 n2 ( \ 2fc + 1 912933 - 9i3'723 / n I n ^ / ^ 

^OQ 933 ^Oq 933 

+ —^773 ^-^-^ —v'^Xoirj) + aDr,Xa{ri) + PvXoiv) + IXoiv) ~ ^iXoiv) = , 

^Oq 933 

where a — a{j, fc), /3 — /3(j, fc), 7 = 7(j, fc) are of the form 

(4.26) a = a(2fc + l)i/2(2j + l), /3 = ^(2fc + l)i/2(2j + 1), 
and 

(4.27) 7 = 71 (2j + 1)2 + 72(2fc + 1)2 + 70, 
where a, /3, 70, 71 and 72 are some explicit constants. 



The proof of Lemma 14.11 is given by a long routine computation and will be 
omitted. 

By Lemma 14. 1[ A4 is an eigenvalue of the second order differential operator 

2fc + 1 922933 - 9^3 n2 2fc + 1 912933 - 913923 / n , n \ 

20o 933 ' 265 933 

2k + 1 qnq33 - qfs 2 , n , « , 
+ ^-73 -ri + aD,, + Pii + J . 

2t)o 933 

This operator is a globally elliptic operator in M. Therefore, it has discrete spectrum 
in L2(M), which is described by the next lemma. 

Lemma 4.2. The eigenvalues of the operator 

n = AD^^ + B{tjD,^ + D„Tj) + Cif + aD,^ +Pv + 1, AC - B^ > , 

are given by 

A„, = (2771 + l)VAC-B^ + -f - -_L_^(Ca2 - 2Ba(3 + Al3^) , 777 e N . 
Proof. The operator H can be written as 



n = XX+ + y/AC - 52 + 7 - \zf , 



where 



^ Mn-r. B + iVAC - 52 

X =A^'^D,, + ^^^ r, + z, 



X^ =^V2 „ , ^^lv^C3^ 



Ai/2 ^ + ""' 
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1 AP- Ba 



"2^1/2 2A1/2VAC-B2' 

and we have 

[X+,X] = 2VAC~B^. 

Therefore, the lemma is immediately proved by following a well-known computation 
of the speetrum of the harmonic oscillator by means of creation and annihilation 
operators. D 



For f)5, we have 












^ _ 2fc + 1 922933 - 9I3 
2foo 933 


,B 


_ 2fc + 1 912933 - 913923 p _ 

262 q^^ 


2fc + 1 91 
26g 


l933 ~ 
933 


9?3 


and a, (3 and 7 are given 


bv (|4.26p and (|4.27p. in particular 


. bv (I2.8|). 


we have 


AC- 


i?2 


(2fc+l)2detQ (2fc + l)2 
4&^ 933 " 462 


2^' 






By Lemma |4.2| the ei 


genvalues of 1)5 have the form 








(4.28) a™ = l^ ( 


2a) 


i/2 

(2fc + l)(2m+l) + j/(j,fc), 


me N, 






y{j.k) 


=7- 


-4(AC-B2)(^«^-2^"^ + ^^^) 








= l'T. 


,{2k + If + uii{2j + if + u^ , 









with 



and 



1^22 = 72, I'D = 70, 



a / 1 911933 - 9?3 -2 , o'?12933 - 913923 _ 5 , , 922933 - 9^3 52 

1^11 = 71 - 3 7 " + 2 ap + oo /3 

a \O0 933 933 933 

Taking A4 as the m-th eigenvalue of ()5 and xo as the corresponding eigenfunction, 
we obtain that the condition (|4.24p holds. Therefore, we can find a solution V2 to 
(|4.23p . Then the condition (|4.22p holds, that allows us to find a solution M4 to (|4.5p . 

Thus, for any j^k,m g N, we have proved the existence of a solution ui,£ — 
0, 1,2,3,4, to the system of equations (|4.ip - (|4.5p with a suitable choice of Xi,£ = 
0, 1,2,3,4 given by (031), BIB . (HIB . B:^ and (IT^ . Setting 



7.'^^ 



€=0 «=o 

we obtain 

The functions it^ have sufficient decay properties. Therefore, by changing back to 
the original coordinates and multiplying by a fixed cut-off function, we obtain the 
desired functions 4''^ i, „^, which satisfy (jl.Sp with 

4 
h sr^ li- 

£=0 
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where 

fJ-j,k,m,e = ^e, ^ = 0,1,2,3,4. 

The functions uq for different j, k and m are orthogonaL Since each change of 
variables, which we use, is unitary, this imphes the condition ()1.4[) . 

5. Periodic case and spectral gaps 

In this section, we apply the previous results to the problem of existence of gaps 
in the spectrum of a periodic magnetic Schrodinger operator. For related results 
on spectral gaps for periodic magnetic Schrodinger operators, see [5] and references 
therein. 

Consider the Schrodinger operator with magnetic potential in the entire Eu- 
clidean space R^ 

H'^ = {hDx, - Ai{X)f + {hDx, - A2{X)f + {hDx, - A3{X)f . 

We assume that the vector magnetic field B = {Bi, i?2, S3) is F-invariant for some 
cocompact lattice F C M^. 
Put 

bo = mm{\B{X)\ : X e R^} 

and assume that there exist a (connected) fundamental domain Q for F and a 
constant eg > such that 

\B{X)\>bo + eo, xedQ. 

We will consider the operator H^ as an unbounded self-adjoint operator in the 
Hilbert space L'^iM?). Using the results of [6], one can immediately derive from 
Theorem 1 1.2 1 the following result on existence of gaps in the spectrum of i/'' in the 
semiclassical limit. 
Theorem 5.1. 
Suppose that 

bo>0, 

and that there exists a unique minimum Xq G Q such that \B{Xo)\ = bo, which is 
non- degenerate: in some neighborhood of Xo 

C~'\X - Xo\^ < \B(X)\ -bo< C\X - Xof . 

Then, for any natural j , k and N, there exist hj^k,N > 0, Cj,^ and Cj^k.N such that 

1 f d^'/^ 



C,,,. > c,, + - ^- j (2fc + l)iV, 

and the spectrum of H^ in the interval [Aj_k.,Bj,k,N\ where 

Aj^k = (2fc + l)bo + (2j + 1)(2A: + l)^l^a^l^h^/^ + h^Cj^u 

and 

Bj.k.N = (2fc + l)bo + (2j + l){2k + l)i/2ai/2/j3/2 _j_ h^^j^^^^^ 

has at least N gaps for any h e]0, hj,k.N] ■ 
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Appendix A. Hermite polynomials 
For a fixed A > 0, put 

Here we gather some well-known formulae, concerning to the Hermite functions /i,„, 
which we need in the paper: 

fhm{t) =j^[h,n+2{t) + {4m + 2)h„,{t) + 4m{m - l)h„^-2{t)] ■ 
^^h^n{t) =^y^[hry^+z{t) + (6m + Q)h„,+i{t) + I2m^h„,^i{t) 



+ 8m(m - l)(m - 2)/i™_3(i)] . 
Dth^{t) ^-\^'^i[h„,+i{t) - 2m/i,„_i(i)] . 

{{tDt + Dtt)h^,h^)^Q. 



{Dtt^ + fDt)Knit) = -^^tiK^+ait) + (2m + 2)h„,+i{t) - 4m^h,n-i{t) 

— 8m(rn — l)(rn — 2)hm-3(t)] . 



{t^Dt+Dtt^)Knit) 

= Y^l^«[^m+5(0 + (6m + 12)h„,+3{t) + 4(2m2 + 4m + 3)/i,„+i(i) 

- 8(2m^ + l)mhrn-i{t) ~ A8m{m - lf{m - 2)/i,„_3(i) 

- 32m(m - l)(m - 2)(m - 3){m - 4)/i,„_5(t)] . 
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